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Emergence of a Kondo singlet state with the Kondo temperature well beyond 1,000K 
in the proton-embedded electron gas: Possible route to high-Tc superconductivity 
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Hydrogen in metals has attracted much attention for a long time from both basic scientific and 
technological points of view. Its electronic state has been investigated in terms of a proton embedded 
in the electron gas mostly by the local density approximation (LDA) to the density functional theory. 
At high electronic densities, it is well described by a bare proton H'*" screened by metallic electrons 
(charge resonance), while at low densities two electrons are localized at the proton site to form a 
closed-shell negative ion H“ protected from surrounding metallic electrons by the Pauli exclusion 
principle. However, no details are known about the transition from H'^ to H“ in the intermediate- 
density region. Here, by accurately determining the ground-state electron distribution n{r) by 
the use of LDA and diffusion Monte Carlo simulations with the total electron number up to 170, 
we obtain a complete picture of the transition, in particular, a sharp transition from short-range 
H+ screening charge resonance to long-range Kondo-like spin-singlet resonance, the emergence of 
which is confirmed by the presence of an anomalous Friedel oscillation characteristic to the Kondo 
singlet state with the Kondo temperature Tk well beyond 1,000K. This study not only reveals 
interesting competition between charge and spin resonances, enriching the century-old paradigm of 
metallic screening to a point charge, but also discovers a long-sought high-Tx system, opening an 
unexpected route to room-temperature superconductivity in a Kondo lattice made of protons. 

PACS numbers: 71.10.Ca,75.20.Hr,71.15.Mb,74.70.Tx 


I. INTRODUCTION 


Physics in the heavy fermion superconductors has been 
understood by the concept of quantum criticality in a sys¬ 
tem of regularly arrayed dense Kondo impurities (Kondo 
lattice)i“— and the spin-fluctuation mechanism is believed 
to be responsible for superconductivity, as inferred from 
the strong correlation between the superconducting tran¬ 
sition temperature and the Kondo temperature T]^—. 
More specifically, Tc is of the order of O.ITk, as shown in 
Fig. [T] plotted based on the information available in the 
literature, from which we can conceive an idea that high- 
Tc superconductivity will be obtained if we can discover 
a Kondo system with very high Tk- In fact, the recently- 
discovered plutonium compounds such as PuCoGas with 
Tc = 18.5K and Tk ~ 260Ki^ may be regarded as the 
successful realization of this idea. Thus we should make 
further pursuit of this idea by searching for a new class 
of Kondo systems with Tk higher than 1,000K. Theo¬ 
retically, this search can be done by the first-principles 
quantitative determination of Tk for the composite sys¬ 
tem of an impurity atom embedded in a metal. 

As first suggested by Debye and HiickeliS, an atomic 
nucleus charge -\-Ze in a metal is screened by accumula¬ 
tion of metallic electrons which is regarded as a charge 
resonance and well described by the linear response the¬ 
ory. This concept of metallic screening is prevailing for 
a century, but because Z is not infinitesimally small, we 
need to consider nonlinear effects in the screening, in¬ 
cluding the contribution from spin fluctuations. The spin 
contribution will be enhanced, if Z is an odd integer, such 
as ^ = 1 (case of a proton), in which a spin-polarized 
bound state might appear at the impurity atom. 
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FIG. 1: (Color online) Superconducting transition tempera¬ 
ture Tc versus Kondo temperature Tk (a characteristic energy 
scale for spin fluctuations) in heavy fermion superconductors. 


With the above basic scientific issues in mind, we have 
concerned with the problem of hydrogen impurity in met¬ 
als which attracts long attention from a technological 
point of viewii, such as hydrogen storage in solids, sen¬ 
sor applications, and catalysis. Its electronic state has 
been investigated in terms of a proton immersed into an 
interacting many-electron system plus a compensating 
background (electron gas: EG) since 1970s, because this 
is an ideal system to study important topics related to 
an impurity in metals, such as the embedding energy, a 
key quantity in the effective-medium theoryi^^— . It is 
also studied from a motivation to improve on the local 
density approximation (LDA) to the density functional 
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theory (DFT) by using the electron distribution n(r) ob¬ 
tained by quantum Monte Carlo methodsi^^—. Nonlin¬ 
ear metallic screening is another extensively examined 
topic in this systenti^^— , but no serious attention has 
been paid so far to the spin resonance effect. 

The homogeneous EG with the average density no 
is specified by a single parameter defined by = 
(3/47rno)^/^ in units of the Bohr radius ob- (We use 
atomic units hereafter.) Its characteristic energy is the 
Fermi energy ep, given by fc|/2 (= 1.84rJ^) with /cf(= 
1.92rJ^) the Fermi momentum (Fig.[5Ka)). Hydrogen, on 
the other hand, has two typical energies (Fig. [2Kb)), the 
Is level eis(=0.5) and the electron affinity eA(=0.0278). 
Ratio of ef to eis or ea determines the character of the 
ground state in the proton-embedded EG; for high no cor¬ 
responding to ^eis (or Ts -C 2), the Is level is buried in 
the continuum of EG and thus no electrons are bound to a 
bare proton H+, leading to a charge resonance (GR) state 
in which H+ is dielectrically screened by accumulation of 
itinerant electrons near the Fermi level /i (Fig.l^Jc)). For 
low no with £f ^ (or » 10), H+ captures two 
antiparallel-spin electrons to form H“. This closed-shell 
negative ion resides in EG with repelling other electrons 
owing to the Pauli exclusion principle (Fig. |2Kd)), but if 
Ef increases and reaches as high as ea , the Fermi pressure 
from EG to the ion becomes so large that the localized 
electrons in H“ begin to spill out into EG. Then a crucial 
question is whether this state at ep ~ £a is the same as 
that in Fig. |2Kc) or not. 


(a) Electron gas (EG) (b) Hydrogen (H) 

Fermi energy: eF=1.84/r,2 ^ 

8,= 0.5 



Eb or 

Is 


l5 



H Affinity: 
£a= 0.0278 

k~ 

1&, +8* 

I IS A 

% _ 

H 


FIG. 2: (Color online) Proton-embedded EG with character¬ 
istic energies in (a) and (b). Three possible ground states 
are schematically illustrated in (c)-(e) corresponding, respec¬ 
tively, to a bare proton H'*" screened by metallic electrons 
with the screening length ~ (a GR state), a closed-shell 
ion H“ confined in EG, and an SSR state in which if the 
Is (virtual) level of H is temporarily occupied by a single 
up- (down-) spin electron, down- (up-) spin clouds of itin¬ 
erant electrons are formed around H for spin screening with 
the screening length kp^. These temporary states are 

superposed with interchanging the roles of spins to make a 
resonance state fl/ssR. 

Intuitively, for ea ^ £f ^£is) we can imagine a spin- 
polarized state made of a single electron with either up 
or down spin captured by H"^, but in view of the con¬ 
cept of spin screening to form a Hondo singletii^ in 


the impurity Anderson model (IAM}2^, we anticipate the 
emergence of not a spin-polarized but a Kondo-like spin- 
singlet resonance (SSR) state 'I'ssR (Fig. HKe)). Because 
there is no clear distinction between conduction and lo¬ 
calized electrons, this SSR state is composed of only 
itinerant electrons near /r without a local spin moment, 
similar to GR, but an important difference exists in the 
screening length; for CR, it is the Thomas-Fermi length 
Atf but for SSR, the Kondo-screening length 

is much longer than fcp leading to an anomalous Friedel 
oscillation2^“— . Then the main aim of this paper is to 
confirm this conjecture about the emergence of SSR in 
the proton-embedded electron gas with determining Tk 
from first principles, but this confirmation is not an easy 
task due to the existence of various difficulties, as we shall 
explain below in some detail. 

The Hamiltonian for lAM, Ha, is written as^ 

Ha — ^ T Hd E ^da^da- + U C^^Cd^C^^^Cdl 

h(T <J 

+ Hcifc(c|;^Cfccr -F cj[,^Cdcr)) (1) 

ha 

in second quantization with use of the annihilation op¬ 
erator Cka for a conduction electron with wave vector k, 
spin cr and one-body band energy Ek, while Cda is an op¬ 
erator to destroy a localized d electron with spin a at the 
energy level Ed ■ The d electrons not only interact to each 
other at the localized site with the strength U but also 
hybridize with the conduction electrons with the strength 
Vdk- Based on Ha, ^'ssr is given as21i2i 

'f'ssR = ^ afe(c^^cfc^ + c];|Cfc;)<i)o, (2) 

k 


where 4>o is the Slater determinant made of conduction- 
electron orbitals and the parameters, ag and a^, are to 
be determined variationally. As Eq. © clearly shows, 
'I'ssR is not descibed by a single Slater determinant but 
is a correlated many-body state. It must also be noted 
that continuum conduction states around the Fermi level 
are indispensable for the construction of this 'I'ssr- 
In our work, we are not allowed to employ this sim¬ 
ple model Ha- Instead, we have to start with the first- 
principles Hamiltonian H, described in atomic units as 


h=-VE 

^ 2 


-E 


Ti-Vi 




Vexti'>’'i)EC^, (3) 


in first quantization. Here we have considered a neutral 
atom of atomic number Z at the origin of coordinates 
immersed into the jellium sphere of radius R and average 
density ng. The number of electrons contained in the 
jellium sphere is 47ri?^no/3= {R/tsY, so that the total 
electron number N is equal to Z + [R/vs)'^, satisfying 
global neutrality, from which we obtain R = {N—ZY^^rs- 
In Eq. the external potential working on an electron 
Vext(f’) is composed of the potential from the nucleus and 
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that from the positive background, written as 


Ve.t{r)=- ^-^- ~ 


N-Z 


e{\r\-R), 


( 4 ) 


with 9(x) the Heaviside function and Cfj represents 
the Coulomb self-energy stemming from both nucleus- 
background and intra-background interactions, given by 

3Z(iV-Z)2/3 ^ 3(iV-Z)5/3 
^^"2 +5 

In solving Eq. ©, we impose the fixed boundary condi¬ 
tion to make the wave function vanish at \ri\=R. Irre¬ 
spective of whether we include the constant term Cfj in 
H or not, there is no problem of divergence in this finite- 
N system, but is needed in order to achieve the mu¬ 
tual global cancellation in energies between the electron- 
background attractive potentials and the repulsive po¬ 
tentials of both electron-electron and intra-background 
interactions for the bulk (N^oo) system22. 

By comparing H in Eq. with Ha, we readily see 
their differences; first, there is no predetermined local¬ 
ized state in H and thus a (virtual) localized state, if 
any, must be determined in the first place in constructing 
'I'sSR with related parameters such as Ed, U, and Vdk, if 
necessary. Second and more importantly, the long-range 
Coulomb interaction, which is completely neglected in 
Ha, works among metallic electrons in H, making solu¬ 
tion of the problem quite difficult. In fact, even in the 
homogeneous EG without the embedded atom, the prob¬ 
lem is sufficiently complicated in both variational^Sr— 
and Green’s-functio n33i34 approaches. Incidentally this 
long-range Coulomb interaction is the source to bring 
about CR and therefore it is indispensable for discussing 
competition between CR and SSR. Since this discussion 
constitutes another important aim of this paper, we can 
never neglect this long-range Coulomb interaction, mak¬ 
ing all theoretical and computational techniques devel¬ 
oped so far for Ha useless to H. 

Usually, the first-principles Hamiltonian is solved by 
either diffusion Monte Caro (DMC) simulations or DFT- 
based methods. The former is an excellent method to 
obtain fairly accurate results for the ground state, but 
it can never directly treat 'I'sSR, because DMC simula¬ 
tions can be done only for finite-fV systems in which all 
levels are discrete, while in constructing 'I'sSR: we need 
continuum conduction states which are allowed only in 
the bulk (iV —>• oo) system. The latter methods can eas¬ 
ily treat the bulk system, but the ground-state physical 
quantities are calculated in terms of a single Slater deter¬ 
minant made of Kohn-Sham (KS) orbitals introduced in 
DFT, so that it is not clear at all as to how much the ob¬ 
tained quantities reflect the highly correlated many-body 
nature of 'I'sSR and how accurate they are, especially be¬ 
cause in actual calculations we always have to resort to 


some approximation to the exchange-correlation energy 
functional E^'^[n{r)] such as EDA. 

Faced with those difficulties, we have decided to focus 
on n{r) rather than the wave function 4 'ssr itself, mainly 
because DFT can, in principle, provide exact n{r) and 
the corresponding ground-state energy Eq by projecting 
the real interacting many-body system to a fictitious aux¬ 
iliary non-interacting system in which n(r) can be calcu¬ 
lated with use of a single Slater determinant made of 
KS orbitals, even if we know nothing about 4 'ssr in the 
real system. This nontrivial assertion, one of the cen¬ 
tral theorems in DFT, is rigorously proved as long as 
the ground state is non-degenerate^, as is the case for 
SSR. Of course, information obtained only through n(r) 
and Eq is limited and useless for discussing transport^ 
and excited-state properties, but we claim that it is still 
plenty enough for our purpose of distinguishing between 
CR and SSR states and determining Tk in SSR. In the 
context of DFT, given exact n{r), it is also an interest¬ 
ing issue to clarify how KS orbitals, which are defined 
in the fictitious system and thus have basically no phys¬ 
ical relevance, behave so as to correctly provide n(r) in 
a strongly-correlated state such as SSR. Thus this clari¬ 
fication constitutes an additional aim of this paper. 

In implementing calculations of n(r), we adopt the 
following strategies: (i) In actual LDA calculations, we 
employ the local spin-density approximation (LSDA)^ 
by choosing a spin-dependent exchange-correlation en¬ 
ergy functional E^^[na-{r)] written in terms of the spin- 
resolved electron distribution n„{r), so that we can 
check a possibility of the spin-polarized (i.e., non-spin- 
singlet) ground state by detecting the difference between 
n-|-(r) and ni{r). (ii) We assess the results in LSDA at 
finite-A^ systems in comparison with those in fixed-node 
DMG^ with taking N up to 170, a much larger size than 
those in previous variational Monte Carlo (VMC) cal- 
culationsi^di. (iii) We obtain n(r) in the bulk system 
by the calculation in LSDA and also by extrapolation of 
DMC data to N^oo. The extrapolated results for n(r) 
in DMC are independent of N and thus we can assume 
that they will be free from any restrictions incurred from 
the fixed-node approximation, the only approximation 
adopted in DMC simulations, because the fixed-node po¬ 
sitions in DMC are prescribed by R (and consequently 
by N for given Vg due to R= {N — Z^^^Vs) in the fixed 
boundary condition, indicating that independence of N 
also suggests independence of the postulated node posi¬ 
tions. (iv) We check whether the obtained n{r) at N —>'00 
exhibits the behavior characteristic to the Kondo SSR 
state or not. More specifically, we look for modulation of 
the Friedel-oscillation period, a very important inherent 
property of the anomalous Friedel oscillation, to which we 
call serious attention for the first time in quantitatively 
determining and consequently Tk from first principles. 

In accordance with those strategies, we have investi¬ 
gated n{r) in finite-A^ systems in both LSDA and DMC 
to find spin-unpolarized ground states with a strong size 
effect and a series of magic numbers (10, 60, 170, • • ■) of 
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N at which convergent results are easily obtained. These 
features can be explained in terms of the emergence of 
SSR with its long-range nature of ^k- Its emergence is 
also signaled in LSDA at iV —^ oo for > 2 by the ap¬ 
pearance of a strange shallow bound KS orbital below the 
conduction band bottom with an unusually long binding 
radius. In spite of the large size effect, for each N, n{r) 
in LSDA agrees very well with that in DMC. 


II. CALCULATION METHODS 

A. LSDA in the finite-A system 

Let us consider the neutral system of a single nucleus 
with atomic number Z embedded in the iV-electron jel- 
lium sphere of radius R. Its Hamiltonian H is given in 
Eq. Q. In LSDA to DFT, the KS equation is written as 


1/2 r-t 


•/Isolated H' ionh 



State d: 

Closed-shell 
H" ion 
confinement 


6 8 10 12 14 

rs 


[-V^/2 -h uP(r)] (l)ia{r) = SiaC^iair), (6) 

where Sia and (jjia are the energy level and the normalized 
wave function for KS orbital i and spin cr, respectively, 
and v^^{r) is the KS potential, determined by 

t;^®(r) = Uext(r)-f J dr' (7) 

where v^^{r;[na]) is derived from A™[no-] through the 
functional derivative as 

v™(r; [ria]) = 6E^^[na-]/Sna{r). (8) 


FIG. 3: (Color online) Ground-state diagram in the proton- 
embedded electron gas in {rs,N~^) space with N the total 
electron number, indicating sharp but size-dependent sequen¬ 
tial transitions among CR, SSR, and closed-shell ion H“. Rig¬ 
orously speaking, SSR is defined only at A —^ oo, but the 
states for finite A directly connected to SSR at A —>■ oo are 
also called SSR. In the bulk system, the screening length 
is predicted to diverge at the CR-SSR boundary or « 1.97 
with the change of Ts, signaling the sharp transition. For 
an isolated H“ ion at A = 2, we have used the exact data 
for n{r)^ to determine = 5.04 by averaging the local 
r's(T’)(= [3/47rn(r)]^'^®) over the weight of n(r) itself. 


With use of the lowest-KS orbitals, ncr{r) is given by 

= (9) 

and n{r) is the sum of nf{r) and n4,(r). The spin density 
n„{r) and consequently N^, with N = should be 

determined by the self-consistent solution of Eqs. dS])-®, 
together with the fixed boundary condition 

</’iCT(r)=0, (10) 


By summarizing the results thus calculated, we have 
obtained a ground-state diagram in (rs,N~^) space, 
shown in Fig. |3l exhibiting sharp sequential transitions 
among CR, SSR, and closed-shell H” ion confinement 
states. Contrary to the previous explanation!^, we claim 
that the very shallow bound KS orbital found in LSDA 
at intermediate densities is not a physical H“ but ap¬ 
pears just to describe the long-range change of n{r) over 
in SSR in the form of an envelope density. We also 
find that hydrogen is most stably embedded in EG in the 
SSR region, especially, optimally firmly at Cg ss 4 with the 
Kondo temperature Tk ~2,100K, indicating our success 
in discovering a long-sought high-Tx system. Thus we 
may expect that superconductivity occurs at a tempera¬ 
ture as high as about O.ITk in a metallic hydrogen alloy 
at ambient pressure in which a macroscopic number of 
protons are regularly embedded in a metal in this den¬ 
sity region to form a Kondo lattice. 

In Sec. II, we explain the calculation methods in both 
LSDA and DMC. In Sec. Ill, we show the calculated re¬ 
sults and in Sec. IV we discuss on the obtained results, 
together with their implications and future directions. 
Finally in Sec. V, we give a summary of this paper. 


at |r| = R = (N — By using those converged 

quantities, we can calculate Eo{N,Z) the ground-state 
energy including the constant term by 


Eq (A, Z) EiCT f dr [uext (r) - (T')]n<^ (r) 

2(7 a 

,u(r)n(r') 


+ : 


dr dr 


\r — r 


- l-E^'^lny+Cf^. (II) 


B. VMC 

With use of the lowest-A KS orbitals thus obtained, 
we can define the Slater determinant $o(t’i,- ■ •, '(’at), with 
which the trial many-body ground-state wave function 
$(ri,- • rjv) for the VMC calculation can be constructed 
in the Slater-Jastrow type as^ 


$(ri,- • •,rAr)=exp[J(ri,- • •, rjv)]4’o(r’i,-'''’ at), (12) 

where the Jastrow function J(t’i,-• •, rjv) contains the 
terms to describe electron-nucleus correlation ui{ri), 
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two-electron correlation U 2 {ri — rj), and three-body 
nucleus-two-electron correlation rj,ri — rj) as 

='^ui{r,) + '^U2{ri-rj) 
i i>j 

+ '^U3{n,rj,ri-rj). (13) 

i>j 

The actual choice of the forms for ui{ri), U 2 {ri—rj), and 
u^iri^Tj^ri—rj) as well as their optimization is done by 
adopting the CHAMP-code package^ as it is. Then the 
expectation value (A) of an operator A is given by 

= ($|A|$)/{$|$). (14) 

By putting A = n(r) = — r^) in Eq. (imi . we 

obtain n(r) in VMC. 

C. DMC 

Starting with the variationally optimized wave func¬ 
tion $ thus determined, we can further improve on the 
ground-state wave function by considering the diffusion 
equation for 4'(r) in the imaginary time r as 

= (15) 

where Er is the reference energy to be adjusted to Eq 
in the course of DMC simulations by removing the r- 
dependence from the asymptotic form of 'I'(t) at r—>- 00 . 
Note that the formal solution to Eq. (HSl) is written as 

T(r) = ^ (16) 

n 

where is the normalized mutually-orthogonal com¬ 

plete set of eigen functions for E[ with the corresponding 
set of eigen energies {E„}. Then, as long as ('I'o|$) ^ 0, 
the asymptotic r-independent wave function T is reduced 
to the true ground-state wave function 4'o, apart from the 
normalization factor. 

It is appropriate to add a comment on the condition 
of (Tol®) 7 ^ 0 here; by invoking the Anderson’s orthog¬ 
onality theorem^, one may argue that ('I'o|4>) vanishes 
in SSR, but this is not correct for the reasons below; (i) 
both and $ include the effect of the impurity atom, 
while the Anderson’s theorem concerns with the relation 
between the wave functions with and without the impu¬ 
rity. (ii) Due to the presence of the Jastrow factor J, 4) 
is not simply given by the single Slater determinant $0 
on which the Anderson’s theorem is proved, (hi) DMC 
simulations are done for finite iV, while the Anderson’s 
theorem becomes valid only at N -^ 00 . 

In order to avoid the notorious fermion sign problem, 
we employ the fixed-node approximation in DMC sim¬ 
ulations. This approximation may bring about unde¬ 
sirable errors in 4', but we try to minimize them by 


seeking for A^-independent results by exploiting the fact 
that the node positions depend on N in the fixed bound¬ 
ary condition, leading to the hope that unphysical node¬ 
position dependent effects will be removed by extracting 
the iV-independent results. In performing actual fixed- 
node DMC simulations at a fixed N, we adopt CHAMP 
again to obtain the stably converged asymptotic wave 
function 4'(ri,- • ^r^r). Then n{r) in DMC is estimated 
by a second-order approximation to the exact expecta¬ 
tion valued, which amounts to 

„(^)DMC ^ 2 (4'|h(r)|$)/(4'|$) - n(r)V“C, (17) 
where n(r)^^^ = (n(r))^^^. 

D. LSDA in the bulk system 

Contrary to VMC and DMC, LSDA allows us to di¬ 
rectly treat the bulk {N —>■ 00 ) system, in which al¬ 
most all states in KS orbitals are continuum ones, for 
which we may write et = fc ^/2 with momentum k and 
= Rki< 7 {r)Yim{r/r) with the spherical harmonics 
Yimirlr) and the radial wave function Rkia-{r) satisfying 
the following boundary condition at r (= |r|) —>-oo : 

Rkicrir) cos[6ia{k)]ji{kr) - sm[dia{k)]ni{kr), (18) 

apart from a normalization factor, where ji{kr) and 
ni{kr) are the spherical Bessel functions and 5ia{k) is 
the phase shift of angular momentum I to be determined 
under the condition of Sio-(oo) = 0 , ensuring that an elec¬ 
tron behaves as a free particle at fc (= |fc|) 00 . There 

is a possibility of finding bound states below the bottom 
of the conduction band (sio- < 0) among KS orbitals for 
which (j>ia{r) —^0 at r —>■ 00 . The Levinson theorem^ 
dictates that the total number of the bound states in 
KS orbitals Nbs is given by J2icri‘^^Yl)Sia{0)/Tr, while 
Z — J2iiTi‘^^Yl)Sia{kF)/TT by the Friedel sum rule^. 


E. Embedding energy 

The embedding energy 6E is defined as the difference 
of the ground-state energies between the atom-embedded 
EG and the system of separated homogeneous EG and 
neutral atom. Usually this concept is relevant only to the 
bulk EG. Thus, denoting the ground-state energy of the 
isolated neutral atom as E^, we can obtain dE through 

SE = lim [Eo{N, Z) - Eq{N - Z, 0)] - Uf. (19) 

N—^OO 

Since SE is of the order 0(1) in comparison with Eq of 
the order 0{N), due care must be exerted in order to 
accurately evaluate SE at N ^ 00 . For this purpose, 
we rewrite Eq. Cl]) with using 5n(r) [= n(r) — no] and 
5[^{k), the derivative of dia{k) with respect to k, into the 
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following form with ensuring the convergence of integrals: 

kp 

la 

j ^ ^ 1 

dr — on(r) + — 
r 2 

“E/ {no])na{r) - <^(r; [no/2])y | 

+ E-^[n^] - E-^[no/2] - E^, (20) 

where BS stands for the set of possible bound states. The 
difference in the contribution of KS energies from contin¬ 
uum states is treated by the consideration of change in 
the density of states, Sl^{k)/'7r. 



dr dr 


f6n(r)dn{r') 


\r — r' 


III. RESULTS FOR PROTON IMMERSION 
A. LSDA for the bulk system 

In line with the previous result^, the ground state 
in the proton-embedded (Z = 1) bulk EG is perfectly 
spin-unpolarized, i.e., n-t'(r) = n^(r) at every r, at least 
for Tg < 15 and our results on n(r) and 6E in LSDA 
are in good agreement with those in previous calcula- 
tionsI^ ~^^d^d^d^ . In Fig. |4l the obtained n{r) normal¬ 
ized by n(0) is plotted as a function of 2kpr/'K for Vg = 1, 
4, and 14, together with the s-wave phase shift Ss{k) 
which is spin-independent and gives by far the largest 
contribution among all 6iaik)s. 

For rg < 1.97, we obtain (5;c,(0) =0 and thus A^bs is zero, 
leading to the typical n(r) in CR with H+ screened by 
metallic electrons in a short range. For Vg > 1.97, on the 
other hand, ^s(O) = tt and thus IVbs = 2 , seemingly im¬ 
plying the sudden appearance of H~ ion confined in EG 
at rs = 1.97— . This must be true, if nL(T’) the localized- 
electron distribution is about the same as n(r) for |r| 
smaller than the H~ ion range, as is the case for Vg > 12.5, 
in which <5s(fc) is distinctive and almost a universal func¬ 
tion of k/k-p, as seen in Fig. SKb). Note that the deep 
dip in n{r) just outside the ion region, as seen for rg = 14 
in Fig. |4l(a), is a typical electron profile describing the 
repulsive action of localized closed-shell electrons to ex¬ 
clude metallic electrons from the ion region by the Pauli 
exclusion principle. 

For 1.97 < Ts < 12.5, however, the “localized” electrons 
behave much differently; at rg = 4, for example, they ex¬ 
tend long up to |r| r:: 82 and concomitantly nL(0) is much 
smaller than n(0), indicating deep and massive penetra¬ 
tion of itinerant electrons into the proton site, but such 
penetration would never be allowed due to the Pauli ex¬ 
clusion principle, if the closed-shell H~ ion were firmly 
constructed. Thus, by remembering that KS orbitals in 
DFT are not necessarily tied with real physical entities 
but just introduced for mathematical convenience to cor¬ 
rectly reproduce n(r), we can assume that this “local¬ 



FIG. 4: (Color online) Electron distribution n(r) normalized 
by n(0) in (a) and spin-independent s-wave phase shift 5s{k) 
in (b) obtained in LSDA for the proton-embedded bulk EG. 
In (a), the dotted curves show the contribution from the KS 
bound state which is absent at rs = 1, must be a real physical 
state representing H“ at Vg = 14, and appears only as mathe¬ 
matical convenience for describing the long-range decrease in 
amplitude of the Friedel oscillation in SSR at Vg = 4. The label 
“localized” with a quotation mark indicates this situation. 


ized” orbital does not represent a real localized state but 
just describes the long-range change of n{r) in the form of 
an envelope density over in SSR. Notice that Sg{kp) in 
this density region lies between 0.7x(7r/2) and 1.3x(7r/2), 
which means that Sg{kp) is close to 7r/2, a value expected 
in the Kondo resonance state in lAM^^. 


B. Comparison between LSDA and DMC 

Before going into a more detailed discussion on SSR, let 
us assess the accuracy of LSDA in comparison with DMG, 
specifically at intermediate densities. At rg = 4, for ex¬ 
ample, in Fig.[^ we see a good agreement between LSDA 
and DMC for n{r) at any r, including the sphere bound¬ 
ary, at each N, though the results in VMC do not match 
so well, assuring the importance to perform DMC for tak¬ 
ing the expectation values in accordance with Eq. dm. 

As for N dependence or the size effect, we find that 
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FIG. 5: (Color online) Examples of the calculated electron 
distribution n{r) obtained by DMC, VMC, and LSDA for (a) 
N = 60 and (b) 170 at = 4. The result in LSDA at A—s-oo 
is also shown. 

N = 170 is not large enough to attain convergence in 
n{r) for |r| >3.5 in both LSDA and DMC. For smaller 
|r|, however, no appreciable difference is seen between 
A = 60 and 170 in DMC (and among N = 60, 170, and oo 
in LSDA) for Cg >4, implying that A = 60 is large enough 
to obtain the convergent n(r) near the proton site. 

C. Cusp theorem and the on-top density 

According to the cusp theorem^, n(r) near a nucleus 
of atomic number Z behaves rigorously in the manner as 

n(r) -ricusp(?') = «(0) exp(-2Zr), (21) 

r ^0 

where ncusp(?') exhibits strictly a linear change with r in 
semilog plots (the dashed lines in Fig. IHKa) and (b)). By 
exploiting this linear behavior, we can rather easily and 
accurately determine the on-top density n(0) in DMC by 
looking at the data in the region of 0.2<r<0.5, in spite 
of the scattered nature of data points for n(r) at r<0.2 
due to the rapid increase in energy scale determined by 
the r“^-Coulomb potential. 

For Vs > 11, n(0) thus obtained is about the same as the 
exact one for an isolated H“ ion^S, as shown in Fig. |6Kb) 
for Ts = 11 — 14, indicating that the state is very close to 
the H~ ion confinement state. A more detailed observa¬ 
tion on the DMC data reveals a systematic change of the 



FIG. 6: (Golor online) Short-range part of n(r) in DMC and 
LSDA at (a) = 4 and (b) Ts = 11 — 14. For comparison, 

the exact result for an isolated negative hydrogen ion H“— is 
also plotted. 


behavior at r <0.2 with the increase of r^; for rg <12.5, 
the majority of data points deviate upward, while oppo¬ 
site is the case for rg>12.5. The upward deviation indi¬ 
cates that the metallic electrons rather easily penetrate 
into the core of the H“ ion, which is not allowed, once 
the closed-shell structure is solidly constructed. On the 
other hand, the downward deviation is consistent with 
the formation of the closed-shell ion, leading to the con¬ 
clusion that the transition to the H~ ion confinement 
state occurs at rgRil2,5, the same rg as that in LSDA. 

In Fig. [3 we plot n(0) by changing rg and A and find 
that the results do not depend on A for rg>3, but they 
do for rg<3, i.e., in the CR-SSR transition region. The 
transition is signaled by a jump in n(0) in both LSDA and 
DMC, although its magnitude decreases with increasing 
A and eventually at A —>■ oo no jump is seen in LSDA 
even at Cg = 1.97 at which the transition is known to 
occur through the abrupt change in i5s(0). 

In an extensive search for favorable A at which we 
can easily obtain the convergent ground state, we find 
that this jump becomes much enhanced at the magic 
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FIG. 7: (Color online) On-top density n(0) plotted as a func¬ 
tion of Ts for A = 60 (and 170 only for Vs = 4) in DMC and 
for N = 60, 170, and oo in LSDA, together with the exact 
results for an isolated H~ ion (A = 2). The transition points 
are indicated by arrows. For the data in DMC, errors are 
within the size of square symbols. 


numbers (10, 60, 170, •••) of A. Let us consider the 
reason for this fact; in LSDA, all KS levels are dis¬ 
crete at finite N and the stacking sequences of the occu¬ 
pied levels are: (ls,2p, 2s), (Is, 2p, 2s, 3d, 4/, 55 ,3p, 3s), 
and (Is, 2p, 2s, 3d, 4/, 3p, 3s, bg, 4d, 6d, 7i, 5/, 8/c, 4p, 4s) 
for N = 10, 60, 170, respectively, reflecting the com¬ 
petition between the —r~^ and potentials in Eq. 0 
and the outermost s orbital, situating very near the Fermi 
level pL at each magic number A, plays a key role in stabi¬ 
lizing the ground state. We recognize that this outermost 
s orbital mimics the SSR state at finite A, at least for 
R < ^K, in view of the fact that the actual SSR state in 
Kondo physics is situated at /i, extending very long over 
the range with an s-wave character. Then the jump in 
n(0) is related to the transition from the empty s-orbital 
state (corresponding to CR) to the s-orbital occupied 
one (corresponding to SSR) with changing A in 8 —>■ 10, 
58 —60 and 168 —170 at rs = 3.85, 2.65, and 2.87, re¬ 
spectively, in LSDA. In DMC, convergent results are also 
easily obtained for the same series of A and the transi¬ 
tion occurs, for example, at Vg ~ 2.65 for A = 58 —>■ 60, 
just as in LSDA. The next A in this series is 340 as¬ 
sociated with the 5s-orbital empty-occupied transition, 
but DMC at this A is currently out of our reach. At 
the second transition into the H“ confinement state, a 
change in the stacking sequence occurs in LSDA. Those 
transition points, along with the second transition point 
obtained in DMC by the data shown in Fig.lHKb), provide 
the ground-state diagram in {rs,N~^) space in Fig. [3] 


D. Embedding energy 

In calculating 6E by use of Eq. (|TO)) at some finite A, 
we need Ao(A —1, 0) in addition to Eq{N, 1), but for an 
even number of A, the ground state of the (A—l)-electron 


system is necessarily a spin-polarized one. Nevertheless, 
for smoothly connecting to the spin-unpolarized ground 
state at A —>• 00 , it is better to calculate in the spin- 
unpolarized situation. In LSDA, by a spin-symmetrized 
sum of KS orbitals by fractionally occupying the or¬ 
bitals at the Fermi level, we can obtain this needed spin- 
unpolarized ground-state energy Ao(A—1, 0) with which 
we can calculate SE for finite-A systems. 
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FIG. 8: (Color online) Embedding energy SE calculated in 
the same situations as those for n(0) in Fig. [3 Errors in 
DMC are within the size of square symbols. 


In DMC, however, we cannot adopt this procedure. 
Thus we first calculate the spin-unpolarized ground-state 
energy Eq{N — 2,0) in DMC, starting with the Slater 
determinant $o(t’i,--- ,rAr_ 2 ) without the outermost s 
orbital. Then we estimate Ao(A—1,0) by 


Ao(A-l,0) = 


A-1 

A-2 


Eo{N-2,0) 


Un- 
(A-l)l/3 


, ( 22 ) 


with the coefficient aAr-i determined through Ao(A—1, 0) 
and Eo{N—2, 0) in LSDA. We have deduced this approx¬ 
imation scheme by considering that the leading term in 
Ao(A—1,0) is, in general, in proportion to (A—1) due to 
extensiveness of the total energy as well as the subleading 
term in proportion to (A—1)^/^ due to the surface-energy 
contribution. In CR in which the outermost s orbital is 
empty even for the system with Z = 1 and A electrons, 
instead of Ao(A —2, 0), we calculate Ao(A, 0) in DMC, 
with which we estimate Ao(A— 1 , 0 ) by a similar strategy. 

In Fig.lHl SE is given as a function of Vg, exhibiting the 
A dependence similar to that in u(0) in Fig. [71 including 
jumps at CR-SSR transitions. Our results for SE agree 
reasonably well with previous onesi^.^— , though due at¬ 
tention to the size dependence was not paid previously. 
These results of SE demonstrate that hydrogen is most 
stably embedded in EG in the form of SSR. Thus the 
concept of SSR is deemed to play a key role in hydrogen 
storage in metals and hydrogen is expected to reside at 
a site with rs ~ 4 in an inhomogeneous metal. 
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E. Anomalous Friedel oscillation 

Basically, the Friedel oscillation due to the presence of 
a proton at the origin is a concept defined in the bulk 
system, but if we try to discuss it with use of the data in 
DMC, we need to eliminate the sphere-boundary effect 
from n(r) obtained in finite-A systems. For this purpose, 
we adopt the following procedure^^; we first calculate 
the charge distribution in the finite-fV system 

corresponding to Eo{N,Z) and then we estimate n{r) 
through the cancellation of the sphere-boundary effect 
by subtracting from it as 


enough information is now available on the node posi¬ 
tions in 5n{r) up to r « 13. Incidentally the outermost 
s orbital contributes much to the oscillation behavior in 
Fig.[^ assuring its importance in the SSR-density region. 

In order to unambiguously confirm the emergence of 
SSR, let us examine the Friedel oscillation in n(r) in the 
light of its general behavior, known as2^“— 


(r)- 


Ano + 




47r2j,3 


cos(2fcF?’ —37r/2-|-2(5^°^)F(r/^K) 


— cos(2fcF?" —37r/2) 


(26) 


n{r) ^ nN,i{r) - +nQ. (23) 

In LSDA, we can employ Eq. (l23ll by using nN-ifi{r) 
which is obtained in the spin-unpolarized situation simul¬ 
taneously with 0), but in DMC, only njv_ 2 ,o(r) 

associated with £’o(IV —2,0) is available. Thus we esti¬ 
mate nj^_ixj(r) from this as 

(24) 

where A is a parameter to correct the difference in the 
sphere radius, defined as A = [(TV —2)/(A—1)]^/^, and 
the r-dependent exponent /3(r) is determined with use of 
the data for the electron distribution in LSDA as 

P{r) = In / In A. (25) 


Here only the s-wave contribution, which indeed domi¬ 
nates others in the present case, is considered and is 
the s-wave phase shift at the Fermi level produced by the 
potential scattering without the Kondo-resonance effect. 
In Eq. (1251) . E(r/^K) = 1 in CR, but in SSR it gradually 
decreases from 1 for r <C to —1 for r due to 

physics of asymptotic freedom^i. Actually in SSR, F{x) 
is known to be a universal scaling function, as explicitly 
given in Fig. ITUI ab 

By appropriately choosing the branch of tan“^ a;, we 
can rewrite Eq. (EH) into 

Air^ 

n{r) -no - , n o cos[2fcFr-|-^(r)], (27) 

with the amplitude A{r) and the phase 6{r), given by 
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FIG. 9: (Color online) Long-range part of n(r) in both 
DMC and LSDA obtained by the elimination of the sphere¬ 
boundary effect. 

An example of n(r) obtained through the above pro¬ 
cedure at Ts =4 is given in Fig. |9l in which we find good 
convergence for r < 7 in changing N by the comparison 
between = 60 and 170, illustrating that this proce¬ 
dure doubles the size-convergent range of r as compared 
with that in Fig. [5] Admittedly, discrepancy is seen in 
the oscillation amplitude between DMC and LSDA, but 
overall good agreement and size convergence are obtained 
in the Friedel-oscillation phase, indicating that accurate 


A{r) = v/l-2F(r/eK) cosi2Si°^) + E{r/, 


(5(r) = tan 


1-F{r/^K) cos(2(5(°^) 
F{r/^K)sin{2Si°^) 


(28) 

(29) 


If is in the range (0,7r/2), 6{r) increases gradually 
from for small r to 5^s'^ E-k/2 for large r in accor¬ 
dance with the change of F{r/^K)- On the other hand, 
if is in the range (7r/2,7r), d{r) decreases gradually 
from for small r to Si^'^—Tr/2 for large r. With taking 
care of such a gradual change in S{r), we can determine 
the ith zero of Sn{r) [=n(r)—no] by 2kpri+Si = iTr+TT/2 
with Si = S{ri) and i = 1, 2, 3, • • •. Then the half pe¬ 
riod of the Friedel oscillation Ai in SSR is given by 
Ai=ri+i-ri = [l-{5i+i -(5i)/7r]AO) with AO)(=7r/2fcF) 
the half period in CR. Because is slightly differ¬ 
ent from Si, Ai is modified from A^) by the amount 
of (^i+i — Si)/TT. This anomalous Friedel oscillation or 
the oscillation-period modulation effect is an important 
consequence of the presence of SSR, but this modulation 
effect has not been well recognized even in the Kondo- 
physics community due to the fact that the effect is to¬ 
tally absent for Si^'^ being equal to a multiple of tt/2 (see 
the denominator in Eq. ((25)) b which happened to be as¬ 
sumed in the previous model calculations^ii^. 

In Fig. fTUI b). we plot Sn(r)r^jZ as a function of 
2kYr j'K at Tg =4 to check whether the modulation effect 
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FIG. 10: (Color online) (a)Universal scaling function charac¬ 
terizing SSR3^*^. (b)Friedel oscillation in SSR in the proton- 
embedded EG at Vs = 4, as seen by the plot of 5n{r)r^jZ 
in both DMC and LSDA to detect the half period of the os¬ 
cillation in comparison with that in CR obtained in the sys¬ 
tem with Z = 0.01, a fictitious tiny charge to ensure that the 
ground state can be accurately obtained in LSDA as a CR 
state even at rs=4. (c)First three modulation factors in the 
half period of the anomalous Friedel oscillation, given through 
Eq. (I29II with use of F{x) in (a), plotted as a function of the 
s-wave potential scattering phase shift <5^°^ for various values 
of kp^K- 


exists or not in our first-principles calculations. In DMC, 
the first four zeros of Sn{r) in units of 2kprijTT, are 
given by 0.994 ± 0.005, 1.970 ± 0.008, 2.949 ± 0.013, and 
3.927 ± 0.019 with the errors estimated by the distribu¬ 
tion of data points around Sn{r) = 0 at each r^. Then 
((5i+i — (5i)/7r the modulation factors for i = 1,2, and 3 
are, respectively, obtained as 0.024±0.005, 0.021 ±0.009, 
and 0.022 ± 0.016 in which the errors are estimated by 
the inclusion of covariance between and r^+i. Those 
results, at least definitely those for i = 1 and 2, assure 
the existence of the shortening of the Friedel-oscillation 


period, confirming the emergence of SSR. 

Once the data for {r^} are known, we can indepen¬ 
dently calculate the modulation factors through Eq. (IMl) 
as a function of and with using F{x) in Fig. lTOT a). 
The results are shown in Fig. llOf c). from which we find 
that our DMC data for the modulation factors (actu¬ 
ally the shortening factors in this case) for f = 1 — 3 
agree very well with those obtained at ~ 35/fcF = 73 
and = 0.86 x (7r/2) (see the horizontal and verti¬ 
cal double-dotted-dashed lines in Fig. [TUT c)). indicating 
that by quantitatively analyzing first few modulation fac¬ 
tors, we can determine both OLnd uniquely, even 
if only the data for the system size much shorter than 
are available. This determination is made possible 
due to the fact that F{r/^K) changes most rapidly at 
r ^k/ 10. By using thus determined and the Fermi 
velocity up, the Kondo temperature Tk is estimated by 
Tk = uf/'Ck ~ 0.0066 Hartree = 2, lOOK with about 10% 
errors, revealing that this is an astonishingly high-TK sys¬ 
tem. With such a high Tk, the system will not exhibit 
the prominent Kondo effects such as the In T anomaly in 
the experiment for T around the room temperature or 
below; it just behaves as a usual Fermi liquid^. 



FIG. 11: (Golor online) Plot of 5n{r)r^ in LSDA in the bulk 
system with both Z = 1 and — 1 at = 4. The contribution 
from the “localized” KS state obtained at Z — 1, 5nj^{r)r^, is 
also given and seen as playing a role of an envelope density. 

The shortening effect in the Friedel-oscillation period 
is also found in LSDA, as seen in Fig[TTl in which 5n{r)r^ 
is plotted for both Z = 1 and —1 in the bulk system at 
Tg = 4. For large enough r outside the SSR binding ra¬ 
dius ^K, the Friedel-oscillation phase for Z = 1 coincides 
with that for Z = —1; the result with Z = —1 is plot¬ 
ted to represent the behavior for the fictitious very tightly 
bound H~ ion confinement state in which the anomalous 
Friedel oscillation is absent. In LSDA, the modulation of 
the Friedel oscillation is brought about by the contribu¬ 
tion from the “localized” density nL(T'). Incidentally, the 
oscillation-period shortening factors in LSDA are much 
smaller than those in DMC; {S 2 — 5i)/tt = 0.0133 and all 
others are less than 0.004, indicating ss 0.014 x {n/2) 
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in LSDA. (We note that Ss^kp) in Fig. |4Kb) is equal to 
(5i°^+7r/2 ss O.SOTtt.) We can easily understand the reason 
for this large difference in between DMC and LSDA; 
because this phase shift is directly connected with the 
wave function, its accurate value will not be obtained 
by LSDA in which the wave function in the fictitious 
non-interacting system is qualitatively different from the 
true correlated SSR wave function. Due to Ri 0 in 
LSDA, we cannot employ the diagram in Fig. [rUf cl to 
determine very accurately. Therefore we estimate 
from the extent of nh{r), which is 82, giving Tk to be 
about 1,900K. In relation to ni,{r), ebs the binding en¬ 
ergy for this “localized” KS state in LSDA is given as 
0.0115 Hartree = 3,600K, which is about twice as large 
as Tk. Thus, although ebs has no direct physical mean¬ 
ing, this quantity seems to be a good measure for the 
magnitude of Tk. In our calculations in LSDA, the val¬ 
ues for ebs are 97 K, 2,200K, 3,500K, 2,200K, and 720K 
for Ts = 2.07, 3, 6, 8, and 10, respectively, indicating that 
in the majority of the SSR region, namely, for 3 < 8, 

we may expect Tk to be well beyond 1,000K. 


IV. DISCUSSION 

Five comments on the present work are in order: 

(i) In SSR, we find a good semi-quantitative agreement 
between LSDA and DMC, but this is by no means fortu¬ 
itous, because this can be understood by the long-range 
nature of ^k which makes the density variation associated 
with the SSR state slow, validating the use of LSDA for 
the calculation of n(r). 

(ii) From our present perspective, we may regard our 
previous study in LSDA on the spin-polarized ground 
states for second-period atoms in Periodic Tabled as a 
successful extension to multi-channel Kondo systems^, 
in which the Hund’s-rule coupling plays a crucial role in 
producing the spin-polarized ground states. 

(iii) In a short term, the immediate next target of re¬ 
search is a hydrogen molecule H 2 immersed in to 

pursue a new concept in chemical bondin g^^i^^ in metals 
in the light of SSR. In fact, Bonev and Ashcroft^ have al¬ 
ready found an interesting bistability between paired and 
unpaired states for Vg > 3.2. In a longer term, we can 
expect fruitful research on new aspects in the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction^ and dilute 
magnetic semiconductors^. 

(iv) If the number of immersed protons is increased up 
to a macroscopic level and those protons are arranged 
into a lattice in a metal with the SSR-density region, 
we may obtain a Kondo lattice with high Tk. Further¬ 
more, if the lattice constant is so arranged as to tune 
the RKKY interaction to be about the same as Tk for 
realizing the quantum-critical situation in the Doniach 
phase diagram^, then as indicated in Fig. [TJ we may 
obtain high-Tc superconductivity in the spin-fluctuation 
mechanism in a metallic hydrogen alloy. In this regard, 


the transition-metal-hydride system such as TiH 2 and 
ZrH 2 ^“— , usually used for batteries, might be a promis¬ 
ing candidate, although it seems that the metallic elec¬ 
tron densities in the metal hydrides so far synthesized are 
too high for our purpose. 

(v) The basic reason for about hundred-times increase 
of Tk and consequently Tc in the hydrogen system com¬ 
pared with those in the /-electron heavy-fermion systems 
is the overall increase of energy scales as seen by the large 
difference in the magnitudes of Is and 4/ energy levels. 
For the same reason of increased energy scales, super¬ 
conductivity with Tc over lOOK has been discussed in 
solid hydrogen at pressures of about 500GPa in the con¬ 
ventional phonon mechanising^— . Very recently, stim¬ 
ulated by the experimental result^, similar discussions 
are made on H 2 S^i^ as well as HsS^ at pressures of 
about 200GPa. Note that the metallic densities in those 
systems are found to be in the CR region, namely, < 2 
(or typically Tg ss 1.4). Therefore, our present proposal 
has nothing to do with those conventional theories for 
superconductivity. We also emphasize that high metal¬ 
lic densities realized at such very high pressures are not 
needed in our proposal, indicating that high-Tc supercon¬ 
ductivity is expected to occur even at ambient pressure, 
contrary to the case of solid H 2 or sulfur hydrides. 

V. SUMMARY 

In conclusion, by employing both DMC and LSDA, 
we have revealed the emergence of SSR in the proton- 
embedded electron gas by the confirmation of an anoma¬ 
lous Friedel oscillation characteristic to the Kondo-like 
spin-singlet state with quantitatively determining Tk 
and emphasized its stability in embedding hydrogen into 
the electron gas. Our work necessitates to modify the 
paradigm of CR in metallic screening to a point charge 
of Z in the sense that, depending on the metallic elec¬ 
tron density and Z, SSR takes the place of CR. This 
work also provides a first concrete example to show how 
the KS orbitals behave to represent the exact n(r) in 
strongly-correlated electron systems such as those de¬ 
scribing Kondo physics, even though they are not al¬ 
ways physically relevant. Finally, this work discloses 
a long-sought high-TK system, opening an unexpected 
route to room-temperature superconductivity in a Kondo 
lattice made of protons and motivating people engaged 
in hydrogen-based physics, chemistry, and technology to 
synthesize hydrogen alloys exhibiting high-Tc supercon¬ 
ductivity at ambient pressure. 
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